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We report an analysis of the interstellar gamma-ray emission from nearby molecular clouds Chamaeleon, R Coro- 
nae Australis (R CrA), and Cepheus and Polaris flare regions with the Fermi Large Area Telescope (LAT). They 
are among the nearest molecular cloud complexes, within ~ 300 pc from the solar system. The gamma-ray emis- 
sion produced by interactions of cosmic-rays (CRs) and interstellar gas in those molecular clouds is useful to 
study the CR densities and distributions of molecular gas close to the solar system. The obtained gamma-ray 
emissivities from 250 MeV to 10 GeV for the three regions are about (6-10) X 10-2^ photons s"! sr"! H- 
atom~^, indicating a variation of the CR density by ~ 20 % even if we consider the systematic uncertainties. 
The molecular mass calibration ratio, Xqq = N(ll2)/WcOt is found to be about (0.6-1.0) X 10^" H2-molecule 
cm"'^ (K km s~^)~^ among the three regions, suggesting a variation of Xco in the vicinity of the solar system. 
Prom the obtained values of XcOt we calculated masses of molecular gas traced by Wco in these molecular 
clouds. In addition, similar amounts of dark gas at the interface between the atomic and molecular gas are 
inferred. 



1. INTRODUCTION 

Observation of high-energy gamma-rays from 
molecular clouds can be used to study the sources 
of cosmic rays (CRs), the CR density, and the dis- 
tribution of the interstellar medium (ISM) in such 
systems. Gamma-rays are produced in the ISM 
by interactions of high-energy CR protons and elec- 
trons with the interstellar gas or the interstellar ra- 
diation field, via nucleon-nucleon collisions, electron 
bremsstrahlung, and inverse Compton (IC) scatter- 
ing. Since the gamma-ray production cross section 
is almost independent of the thermodynamic state of 
the ISM, gamma-rays have been recognized as a pow- 
erful probe of the distribution of interstellar matter 
(e.g., Stecker 1970, Lebrun et al. 1983 and Blemen 
et al. 1984) If the gas column densities are estimated 
with good accuracy by observations in other wave- 
bands such as radio and infrared, the CR spectrum 
and density distributions can be examined as well. 
Molecular clouds that are within 1 kpc from the so- 
lar system (namely nearby molecular clouds) and have 
masses greater than a few lO'^ Mq are well suited for 
an analysis of their gamma-ray emission to investigate 
the distribution of CR densities and interstellar gas 
since they are observed at high latitudes and therefore 
largely free from confusion with the strong emission 
from the Galactic plane. Study of such nearby molec- 
ular clouds in gamma-rays can be dated back to the 
COS-B era (Bloemen et al. 1984) and was advanced 
by the EGRET on board Compton Gamma-Ray Ob- 
servatory (e.g.. Hunter et al. 1994). Although some 
important information has been obtained on proper- 
ties of CRs and the ISM by these early observations, 
detailed studies have only been performed on giant 
molecular clouds with masses greater than '^10^ Mq 
such as the Orion complex (e.g., Digel et al. 1999). 



The Large Area Telescope (LAT) (Atwood et al. 2009) 
on board the Fermi Gamma-ray Space Telescope has 
improved the situation significantly thanks to its un- 
precedented sensitivity. In addition, recent develop- 
ments in studies of the ISM (e.g., Kalberla et al. 2005, 
Dame et al. 2001 and Grenier et al. 2005) allow us 
to investigate the CR spectra and density distribution 
with better accuracy. 

Here, we report an analysis of Fermi LAT ob- 
servations of the Chamaeleon, R Coronae Australis 
(R CrA), and Cepheus and Polaris flare molecular 
clouds. They are among the nearest 300 pc from 
the solar system) molecular clouds. No detailed study 
of CR and matter distributions for the Chamaeleon 
and R CrA regions had been performed before since 
they have rather small masses {'^ 10'' M©) and con- 
sequently small gamma-ray fluxes. We also analyzed 
in detail the region of the Cepheus and Polaris flares 
which was included in the Fermi LAT study of the 
second Galactic quadrant (Abdo et al. 2010a). De- 
tails of this proceedings are presented in a published 
paper (Ackermann et al. 2012a). 



2. DATA ANALYSIS 



For this analysis, we have accumulated events ob- 
tained from 2008 August 4 to 2010 May 9. During 
this time interval the LAT was operated in sky survey 
mode nearly all of the time and scanned the gamma- 
ray sky with relatively uniform exposure over time 
(within 10% in regions studied). We used the stan- 
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dard LAT analysis software, Science Tools^ version 
v9rl6p0 and the response function P6_V3_DIFFUSE, 
which was developed to account for the detection in- 
efBciencies due to pile-up and accidental coincidence 
of events (Rando et al. 2009). We set the lower en- 
ergy limit at 250 MeV to utilize good angular reso- 
lution (68% containment radius is ;< 1.5° above 250 
MeV) and the upper energy limit at 10 GeV because 
of limited photon statistics. The gamma-ray count 
maps obtained {E > 250 MeV) in the three regions 
are shown in Figure [TJ Point sources with high sig- 
nificance (test statistic, TS^, greater than 50) are also 
plotted. 

We modeled the gamma-ray emission as a linear 
combination of maps tracing the column density of the 
interstellar gas. This approach has been successfully 
applied in recent studies of diffuse gamma-rays by the 
LAT (e.g., Abdo et al. 2010a and Ackermann et al. 
2011). With the usual assumptions of optical thinness 
and that CRs uniformly thread the ISM, gamma-ray 
intensity (s~^ cm~^ sr^^ MeV~^) at a given energy 
can be modeled as 



2 

I^{l,b) - 5l9HM-A^(Hi)(Z,6),-t-(7co-l^co(^,^') 
+ qAv ■ Av,.csil,b) + cic ■ Iic{l,b) 

+ /iso+^PSj 

j 

where i labels the regions along the line of sight sepa- 
rated in the analysis (1 for the Gould Belt and 2 for the 
rest), gni.i (s^^ sr^^ MeV~^) and gco (s~^ cm~^ sr~^ 
MeV~^ (Kkms^^)^^) are the emissivity per H i atom 
(traced by the 21 cm line of atomic hydrogen) and per 
Wco unit (as a tracer of molecular gas), respectively. 
9Avres {s~^ cm~^ sr""^ MeV~^ mag~^) is the emissiv- 
ity per Avros magnitude which accounts for the gas 
not traced well by H i and CO: we constructed visual 
extinction (Av) maps providing an estimate of the to- 
tal column densities on the assumption of a constant 
gas-to-dust ratio, and obtained map by the fitting a 
linear combination of the iV(HI) and Wco niaps. lie 
and /iso are the IC model and isotropic background 
intensities (s~^ cm~^ sr~^ MeV""'^), respectively. We 
used GALPROP^ (e.g., Strong and Moskalenko 1998) 
for the IC model, and a publicly available isotropic 



spectrum^ for the isotropic component. PSj repre- 
sents contributions of point sources included in the 
IFGL catalog (Abdo et al. 2010b). See Ackermann 
et al. (2012a) for details of the analysis. 



3. RESULTS 

Figure[2](left) shows the obtained gamma-ray emis- 
sivity spectra per H-atom of each region assuming 
the spin temperature Tg = 125 K. The emissivities 
above 250 MeV are (5.9 ± O.l^tat t?:osys) x lO'^^ 
photons s"^ sr^^ H-atom^^, (10.2 ± 0.4stat ti'ysys) x 
10-2"^ photons s^i sr^i H-atom" \ and (9.1 ± O.Sgtat 
to'esys) X 10"^^ photons s"^ sr"^ H-atom"^ for the 
Chamaeleon, R CrA, and Cepheus and Polaris flare 
regions, respectively. In order to examine the system- 
atic uncertainty due to the optical depth correction, 
we also tried to fit the data with maps obtained by as- 
suming Tg = 100 K and under the approximation that 
the gas is optically-thin. We evaluated the uncertainty 
of the isotropic component to be ±10% by comparing 
the model we adopted and those derived in other LAT 
studies of mid-latitude regions (Abdo et al. 2009a and 
Abdo et al. 2009b). We thus reran the analysis as- 
suming a 10% higher and lower intensity for the fixed 
isotropic component. We also investigated the effect 
on the systematic uncertainty due to the IC compo- 
nent by using different IC model maps which are con- 
structed under different assumptions about the distri- 
bution of CR sources and intensity of the interstellar 
radiation field (e.g., Ackermann et al. 2012b). The 
effects of the uncertainty of the Tg, isotropic compo- 
nent, and IC models are quite comparable, therefore 
we added them. The resultant peak-to-peak uncer- 
tainty of the local H i emissivity is less than ~ 20 % 
across the energy range for three regions investigated. 



4. DISCUSSION 

4.1. CR Density and Spectrum Close to 
the Solar System 

We then compared the obtained H i emissivity spec- 
tra among the three regions, and model spectra^ for 
the local interstellar spectrum (LIS) as shown in Fig- 
ure[2](left). Whereas the spectral shapes for the three 



^Available from the Fermi Science Support Center 
i jhttp: / /fermi. gsfc.nasa.gov/ssc/ ). 

^TS is defined as TS = 2(lnL — lnLo)i where L and Lq are 
the maximum likelihoods obtained with and without the source 
included in the model fitting, respectively; see (Mattox et al. 
1996) 

^http: / /www. mpe.mpg.de/~aws/propagate. html 



^isotropic jem_v02.txt from http://fermi.gsfc.nasa.gov/ssc/ 
data/access /lat/BackgroundModels.html 

^The model is calculated from the LIS compatible with the 
CR proton spectrum measured by Alcaraz et al. (2000) and 
Sanuki et al. (2000), under the assumption that the nuclear 
enhancement factors are 1.45 and 1.84 (Mori 2009); see Abdo 
et al. (2009b). 
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Figure 1: Gamma-ray count maps above 250 MeV for the Chamaeleon (top left), R CrA (top right), and Cepheus and 
Polaris flare (bottom left) regions, smoothed with a Gaussian of a = 0.5° for display. The contours indicate intensities 
Wco of the 2.6 mm line of CO (with the levels of 4, 8, 12, and 16 K km s"^) by Dame et al. (2001), as a standard 
tracer of the molecular gas. White circles show the positions of point sources with high significance (TS > 50) in the 
First Fermi LAT catalog (IFGL) by Abdo et al. (2010b). The inner squares indicate the ROI analyzed for each of the 
regions. Point sources outside of this ROI but inside the outer square are taken into account in the analysis. 



regions studied here agree well with the LIS model, 
the absolute emissivities differ among the three re- 
gions. We note that the systematic uncertainty of the 
LAT effective area (5% at 100 MeV and 20% at 10 
GeV; (Rando et al. 2009) does not affect the rela- 
tive value of emissivities among these regions. The 
effect of unresolved point sources is small, since we 
have verified that the obtained emissivities are robust 
against the lower threshold for point sources between 
TS — 50 and 100. Although the emissivities of the 
R CrA region and the Cepheus and Polaris flare re- 
gion are comparable, that of the Chamaeleon region 
is lower by ~ 20%, even if we take the systematic 
uncertainties into account. Therefore, the LAT data 
point to variations in the CR densities within 300 
pc in the Gould Belt. If the CR density has a varia- 
tion by a factor of ~ 1.2 in the neighborhood of the 
solar system, this requires a serious reconsideration 
of a smooth CR density often adopted for simplicity, 
and may have an impact on the study of the CR source 
distribution and diffuse gamma-ray emission. We note 
that CR sources are stochastically distributed in space 
and time, and this may produce a CR anisotropy de- 



pending on the propagation conditions (e.g., Blasi & 
Amado 2012a and Blasi & Amato 2012b). Study of 
other regions and more detailed theoretical calcula- 
tions are needed for further investigating this issue. 
On the other hand, we need to further investigate 
possible systematic uncertainties due to variations in 
-'^co (see the next paragraph) or dust-to-gas ratio we 
assumed in the construction of the Avj-cs map. 

4.2. Molecular Masses in the Interstellar 
Clouds Studied 

Under the hypothesis that the same CR flux pene- 
trates the H I and CO phases of an interstellar com- 
plex, we can calculate the molecular mass calibration 
ratio, Xco, as Xco = Qco/(2qhi)- As shown in Fig- 
ure [5] (right), we derived the Xco values from a lin- 
ear fit of the (qhi, Qco) points. The obtained Xqo 
values are summarized in Table |T| (0.96 ± 0.06stat 
t°:i^,ys) X 1020 cm-2 (K km s-^)-\ (0.99 ± 0.08,tat 
t^ofosys) X 1020 cm-2 (K km s-i)-i, and (0.63 ± 
0.02stat to.o?sys) X 10^0 cm-^ (K km s'^)-^ for the 
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Figure 2: (left): Hi emissivity spectra of each region compared with the model for LIS with the nuclear enhancement 
factors of 1.45 and 1.84 (Mori 2009). The shaded areas indicate the systematic uncertainty: vertical lines 
(Chamaeleon), diagonal lines (R CrA) and the gray area (Cepheus and Polaris flare), (right): CO versus Hi 
emissivities. Each point corresponds to an energy bin. Errors are statistical only. 



Chamaeleon, R CrA, and Cepheus and Polaris flare 
regions, respectively. The obtained value of Xqq for 
the Cepheus and Polaris flare region is ^ 20 % lower 
than that reported by Abdo et al. (2010a). Abdo et 
al. (2010a) includes in their study also the Cassiopeia 
molecular cloud in the Gould Belt, and due to different 
ROIs considered, the qhi emissivity was also different. 
Xqo of the Chamaeleon region is similar to that of 
the R CrA region, whereas that of the Cepheus and 
Polaris flare region is ^2/?, of the others. The LAT 
data thus suggest a variation of Xqq on a ~ 300 pc 
scale. 

Using these Xqq values, the molecular masses 
traced by Wco can be calculated. The mass of the 
gas traced by Wco is expressed as 

where d is the distance to the cloud, mu is the mass 
of the hydrogen atom and /i = 1.36 is the mean 
atomic weight per H-atom (Allen 1973). From this 
equation the mass of gas traced by CO is expressed 
as Mco in Table HI we obtained 5x10^ Mg, ~ 
10^ Mq, and ~ 3.3x10"* Mq for the Chamaeleon, 
R CrA, and Cepheus and Polaris flare regions, respec- 
tively. From the linear relation between qui and qawcs 

vrcs vrcs 

• qm) we can calculate the masses 
of additional gas traced by gAvrcs with a procedure 
similar to that for CO and results are summarized 
in Table HI We thus obtained mass estimates for the 
Chamaeleon and R CrA regions similar to previous 
ones (Dame et al. 1987) if we consider the total mass 
(traced by Wqo and Avrcs), valthough the procedure 
is not straightforward since the gas traced by Avj^os 
is extended in a much larger region of the sky. De- 
tailed study of the matter distribution in the interstel- 
lar space by comparing gamma-rays and other tracers 
will be reported elsewhere. 



5. Summary and Conclusions 



We have studied the gamma-ray emission from the 
Chamaeleon, R CrA, and Cepheus and Polaris flare 
molecular clouds close to the solar system (< 300 pc) 
using the first 21 months of Fermi LAT data. Thanks 
to the excellent performance of the LAT, we have ob- 
tained unprecedentedly high-quality emissivity spec- 
tra of the atomic and molecular gas in these regions 
in the 250 MeV - 10 GeV range. 

The gamma-ray emissivity spectral shapes in three 
regions agree well with the model for the LIS (a model 
based on local CR measurement), thus indicating a 
similar spectral distribution of CRs in these regions. 
The emissivities, however, indicate a variation of the 
CR density of 20 % within ~ 300 pc around the 
solar system, even if we consider the systematic un- 
certainties. We consider possible origins of the vari- 
ation are non-uniform supernova rate and anisotropy 
of CRs depending on the propagation conditions. 

The molecular mass calibration ratio Xco f^or the 
Chamaeleon cloud and the R CrA cloud are compa- 
rable, whereas that of the Cepheus and Polaris flare 
region is ~ 2/3 of the others, suggesting a variation 
of Xco in the vicinity of the solar system. From the 
obtained values of Xco, the masses of gas traced by 
Wco in the Chamaeleon, R CrA, and Cepheus and 
Polaris flare regions are estimated to be ~ 5 x 10^ 
Mq, ~ 10^ Mq, and ~ 3.3 x 10^ Mq respectively. 
Similar amounts of gas are inferred to be in the phase 
not well traced by the H i or CO lines. Accumulation 
of more gamma-ray data, particularly at high ener- 
gies, and progress in ISM studies, will reveal the CR 
and matter distribution in greater detail. 
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Table I Xco, ^Avrcs and masses in the interstellar clouds in each region 





^co 


Mco 


-^Avrcs 




(X 


10^0 cm-2 (K km s~^)-^) 




(x 10^^ cm"^ mag"^) 




Chamaeleon 


0.96±0.06stat±0.13,y, 


~5xlO^ 


0.22±0.01,tat±0.08,y, 


-2.0x10" 


R CrA 


0.99±0.08,tat±0.14,y, 




0.21±0.01stat±0.02,y, 


~103 


Cepheus & Polaris 


0.63±0.02,tat±0.08sys 


~3.3xl0'' 


0.14±0.01,tat±0.03eys 


-LSxlO"* 
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